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Abstract. We provide an explicit formula for the Tornheim double series 
T(a, 0, c) in terms of an integral involving the Hurwitz zeta function. For 
integer values of the parameters, a = m, c = n, we show that in the most 
interesting case of even weight N := m + n the Tornheim sum T(m, 0, n) can 
be expressed in terms of zeta values and the family of integrals 

\ogT{q)B k {q) Cl, +1 (27rg) dq, 

la 

with k + I = N, where B k (q) is a Bernoulli polynomial and Cl;_)_i(x) is a 
Clausen function. 



/' 

Jo 



1. Introduction 



The function 



oo oo 1 

T (a, b, c):=J2Y, a 67 I W 
z — i t — 4 r a s" (r + s)° 

r—l s—1 v ' 

was introduced by Tornheim in [11 . For a, b, c G R, the series is convergent if 

(1.2) a + c> 1, 6 + c> 1, and a + b + c> 2. 

In the case (a, 6, c) = (to, /c, n), with to, fe, n C N U {0}, we define the weight of 
the Tornheim sum T(to, fc, n) as the positive integer N — to + k + n. 

We have previously derived [7] an analytic expression for the general Tornheim 
sum T(a,6, c) in terms of integrals involving the Hurwitz zeta function £(z,q), 
defined as the meromorphic extension to the whole complex z-plane of the series 

OO -j 

(1.3) C(*,9) 



which is defined for Rez > 1 and q ^ 0, —1, —2, ■ • ■ . Our expression for T(a, 6, c), 
recalled later in Theorem 13. 2[ is valid for a, b, c e R — N U {0}, provided the con- 
vergence conditions (|1.2p are satisfied. Using this result and a limiting procedure, 
we derived similar formulas for T(m, k, n), with to, k, n G N. 

In this paper we derive a formula for the Tornheim sum T(a, 0, c), valid for 
a, c G R — N with a > 2 and c > 2. A limiting procedure will then provide an 
analytic expression for the sums T(m, 0, n), with to, n G N — {1}. 
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These expressions for the sums T(m, 0, n) are of practical interest. In fact, Huard 
ct al [10] established the relation 
(1.4) 

m + k — i — 1\ „,. _ , T i m + k — i — 1 



T(m,fc,n)=E( )T(i,0,N-i) + J2( \T(i,0,N-i), 

i=i v m * ' i=l \ 1 J 

for the Tornheim sum of weight iV = m + k + n. Therefore, it suffices to consider 
sums of the form T(m, 0,n). The convergence of the series requires n > 1 and 
m + n > 2, thus the sum T(m, 0, 1) diverges. 

Introduce the spaces 

(1.5) Z N := {T(m, fc, n) : to, fe, n G N U {0} with n + to > 2, fe + n > 2 

and iV = m + fc + n>3}, 

and 

(1.6) Z° N :={T(m,0,n)eZ N }. 
The following result is contained in formula (| 1 .4[) : 

Proposition 1.1. Every sum in Zjq is a linear combination of terms in Z^ with 
coefficients in N. 

Note 1.2. The sums T(m, 0, n) appearing in (|1.4j) can be written as 

r— 1 s — 1 ^ ' 



E- 



i-i >r 2 



Therefore T(m, 0, n) is a special case of the multiple zeta value (= MZV) 



fe 

-Si 



(1.8) C(si,s 2 ,...,Sfe)= E II 

ri >r 2 > ■ ■ ■ >r k >0 j'= 1 

nameljQ, 

(1.9) T(m,0,n) = C(n,m). 

There is a vast literature on MZV and the reader if referred to Chapter 3 of [2] for 
an introduction to this topic. 

Note 1.3. In the case of odd weight, [10] gives the relation 

(1.10) r(m,0,n) = (-ir E ( m+ l~ 2 l~ 1 )am(rri + n-2j) 

j=0 \ 712 J 

+ (-l) m E f m + " f ~ ^ C(2i)C(m + n 2j) - ±C(m + «)> 



^The double zeta function appearing here should not be confused with the Hurwitz zeta func- 
tion in JT3}. 
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valid for m > 1. Here 

(i-ii) c(*) = E- 

71=1 

is the classical Riemann zeta function. This function has an analytic extension to 
C — {1}, the point ,s = 1 being a simple pole. Recalling that £(0) = —1/2, the 
previous result can also be restated in the following terms: 

Proposition 1.4. Assume the weight N is odd. Then the sums in Z^ can be eval- 
uated as linear combination of the products ((2j)(^(N — 2j), j = max { |_-?rj , L^T^J }j 
with integer coefficients. 

The idea of [10] is to produce a linear system of equations for the unknowns 
Ti := T(i,0,N-i), for 1 < i < N - 2. This system has full rank in the case N 
odd and its solution yields (|1.10|) . Their methods also produce analytic expressions 
for Tornheim sums of small even weight, but they fail in general. In particular, for 
weight N > 8, the system of equation mentioned above is not of full rank. Section 
[6] contains a discussion of this issue. 

For example, the class Z\ contains the sums 

(1.12) r(l,0,3) = ~C(4) and T(2, 0, 2) = -£(4), 
and in the space Z§ we find the four sums 

T(l,0,5) = -ic(3) 2 + ^C(6), 

T(2,0,4) = C 2 (3)-^C(6), 

T(3,0,3) = ic(3) 2 -ic(6), 

OK 

(1.13) T(4,0,2) = -C 2 (3) + -C(6). 
Huard et al. [TU] also gave the relation 

49 

(1.14) 5T(2, 0, 6) + 2T(3, 0, 5) = 10C(3)C(5) - — C(8), 

for the case of weight 8, but they are unable to evaluate the individual terms 
T(2, 0, 6) and T(3, 0, 5). The question of an analytic expression these sums remains 
open. 

In this paper we give particular consideration to the Tornheim sums T(m, 0, n) 
of arbitrary even weight N — m + n. For each even N, only two of these sums 
have known closed expressions in terms of zeta values, namely T(1,0, iV — 1) and 
T(N/2,Q,N/2), 

Tornheim established the result 

(1.15) T(0,0,N) = ((N-1)-((N), 7V>3, 

which appears as Theorem 5, page 308 of Tornheim [11 , and the companion formula 



(1.16) T(l,0,7V-l) = i 



N-2 

(N-l)((N)-yC(i)C(N-i) 



i=2 



N > 3, 



where N = n + 1 , can also be found in [TT] . 
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The main result of this paper is an analytic expression for the Tornheim sums 
T(m, 0, n) of even weight, with m, n > 2, in terms of a family of integrals involving 
the log-gamma function logr(g), the Bernoulli polynomials Bk(q), given by the 
generating function 



~fcT' 

k=Q 



and the Clausen functions CI; (a;), defined as 

OO . T 

(1.18) Cl 2n (z):=£™ neN, 

fc=l 
OO , 

(1.19) CWH,):^^, neNU{0}. 

k=\ 

For example, we obtain 

(1-20) T(6, 0, 2) = IC(8) - 6C(3)C(5) - Y 2 \, 



with 



and where 



Y* 6 := -7T 6 (X 0|6 - 2*i )5 + X 2A ) - 6C(7) log 2^, 



(1.21) X Kl := (-l)LV2j JL [ i ogT ( q )B k (q) Cl i+1 (27rg) dq. 

(27T)' Jo 

In the general case, we show that all the Tornheim sums of even weight N can 
be expressed in terms of values of the Riemann zeta function and integrals of the 
form 



n22 , Y * 2(2^ " w( m \X 

(1.22) Y mtN _ m - m , (iv _ m _ 2)! ^(-l) [jjX^-j 

+ -IT- 1 CiV-llog27r, 
\m — 1/ 

where iV is the weight and m is even in the range 2 < m < 2 [^g^J • 

The rest of the paper is organized as follows. Section [2] contains all the main 
theorems we will prove in the subsequent sections. Section [3] derives the expression 
for the Tornheim sum T(a, 0, c) in terms of an integral involving the Hurwitz zeta 
function, starting from a more general result derived in [7j . Section 3] computes 
the limit of T(a,0,c) as a — > m and c — > n, with m,n S N — {1}. In Section 
[5] we consider the particular case of Tornheim sums of even weight N and show 
that they can be expressed in terms of zeta values and the family of integrals Xk,i, 
defined above, with k + l — N. The explicit evaluation of this last family of definite 
integrals remains a challenging problem. Finally, in Section [5] we give a systematic 
list of evaluations for small even weight. 
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2. Main results 

We start by introducing some auxiliary special functions that play a role in our 
derivations. First, we have the Bernoulli functions 

(2.1) A k (q)=kC(l-k,q), keN, 
introduced in [3 [6] , and the kernel 

(2.2) K (q) = logsuiTrg. 
The first result is established in Section [3J 

Theorem 2.1. Let a.ceR-N with a, c> 2. Then 



(2.3) T(a, 0, c) = 4A(o)A(c) sin (—) 



sin (™) C (l - «)C(1 - c) - C(l-* -c)g( 0; c) 
V 2 J | sv AV ' l-tan(^)tan(f ) j 

\ cos (t) / [c(1 _ a " q) c(1 _ a ' 1 _ 9)1 c(1 _ c ' 9) 



cot 7rg (ig 



where 



(2.4) AW-™ " : 



(2-k) 1 ~ z (2ny- z T(z) sinTrz 

The explicit expression (|2.3p for T(a, 0, c) allows us to consider the limit a — » 
m, c — ► n, for to, n G N — {1}. The value of T(m, 0, n) is found to be given in terms 
of the basic integrals: 

(2.5) I BB (k,l):= [ B k {q)B l (q)K(q)dq, 



1 



I AB {k,l):=- A k (q)B l {q)K(q)dq, 

T JO 
1 



I AA (k,l):=— A k {q)A l {q)K{q)dq, 
n Jo 

1 Z -1 

Jaa(M) := — / ^(9)^(1 - q)K{q) dq. 
1X " Jo 

where B k (q) is a Bernoulli polynomial and A k (q) and i^g) are the functions intro- 
duced in (|2.ip and (|2.2p . respectively. 

Theorem 2.2. Let to > 2, n > 2 e N. The Tornhcim sum T(m, 0, n) is given by 

(2.6) T(m, 0, n) = C(m)C(») - |C(m + n) + (-1) 1 ; , , l a (m, n), 

m\n\ 

where: 

(a) to and n even: 

(2.7) £2(171, n) = mlAB(m — l,n) + nlAsim, n - 1), 
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(b) m and n odd: 

(2.8) liim, n) = mlAB(n, m — 1) + nlAB(n — 1, m). 

(c) m odd and n even: 

(2.9) £ 2 (m,n) = ^(mI BB (m- 1, n) + nI BB (m, n - 1)) , 

(d) m even and n odd: 

(2.10) £2(171, n) = —ml AA(m ~~ lj n ) — nlAAim, n — 1) 

— mJAAi'm — l,n) + tiJaa ( m , i — 1), 
The proof of this theorem is the subject of Section 0J 

In the case of even weight m + ji = N, the Tornheim sum T(m, 0, n) is given 
in terms of the single family of integrals lAB{k,l). The next theorem shows that 
these can be expressed in terms of the family of integrals 

(2.11) X ktl := (-l)L'/aJ " f\ogT( q )B k (q)Gh +1 (27Tq)dq. 

Theorem 2.3. Assume m, n G N satisfy m, n >2 and that the weight JV := m + n 
is even. Dehne 

(2.12) r!(m,n):=C(m)C(n)-icW 
and 

(2.13) T 2 (m,n) : 
where 



N/2-2 



£ 2 _ 2fc )c(2fc+i)C(^-i-2fc)+(-ir/ 2 -X* i , 

fe=i v TO / 



2(27r) jv ~ 2 m - /m\ 



m!(iV - m- 2)! ^ v y \ j 
( 2 - 14 ) + (-l)f- 1 (^: 2 ) C (^-l)log2. 
F* -F^ + C-^T- 1 C(iV-l)log27r. 



771—1 



Then T(m, 0, n) can be written as 
(2.15) T(m,0,rc) 



Ti(m, n) + Ti(m, n), m and n odd, 
Ti(to, n) + T2(n, m), m and n even. 



It follows that, for even weight N — m + n, T(m, 0, n) is determined by either 
Ym,n or y„ jm , depending on the parities of m and n. Note that both Y m ^ n and Y n<m 
are linear combinations of integrals X^j with fixed k + I = N — 2 = m + n — 2. 

We discuss this theorem in Section [5l 
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3. An expression for the Tornheim series T(a, 0, c) 

In this section we present an analytic expression for the Tornheim double series 
T(a, 0, c), valid for a > 2, c > 2 and a, c g N. Note that T(a, 0, c) will be finite if 
a > 1 and c > 1. We employ here and in Section [4] the shorthand notation 

(3.1) C(z,q) :=C(1 -*,<?) and CO*) == C(l - *). for z ^ 0. 

All the properties of the Hurwitz zeta function that we will use can be expressed 
in terms of the function ((z, q) as follows: 

(3.2) ±az,q) = (z-l)az-l,q), 

aq 

(3.3) t(k,q) = ~B k (q), 

(3.4) Z'{k,q) = ~A k {q). 

where Bk(q) is a Bernoulli polynomial and A k (q) is the Bernoulli function (|2.ip . 

The identities derived below appear from integration by parts. The restrictions 
imposed on the parameters guaranteee that the boundary terms vanish. Recall that 
the Hurwitz zeta function satisfies the identity 

C(z,q) =q- z + C{z 1 l + q) 1 

which implies Q{z, 0) = C{z, 1) if z < 0. Equivalently, 

(3.5) C(*,0)=CM) if z>l. 
Theorem 3.1. Let a, c e M - N with a, c> 2. Then 

(3.6) T(a,0,c) = 4A(a)A(c)sin 

(f){c»c»- T « a + c ^ a ^ 1 



sm i 



tan (f) tan (f) j 
C ° S (t) / ^ a,qS) ~ ^ ai1 ~ ^ C{c,q) cot irqdq 



1 /"ira 
~ 2 

The proof is obtained by analyzing the behavior as e — > of T(a, b, c) given in 
below. This was first derived in [TJ. The parameter b is changed to e in order 
to remind ourselves that it is small. 

Theorem 3.2. Let a, b, c S R, satisfying a + c > 1, 6 + c > 1, and a + b + c > 2, 
and define the auxiliary function A(z) as in (|2.4[) . 

fttf AM •= r ( 1-z ) - g 

k ; lJ ' {2iry- z (2 7 r) 1 - z r(z) sinTrz' 

For o, d,c^ffwe have 

(3.8) T(a, b, c) = 4A(a)A(6)A(c) sin (y ) Q(a, b, c) 

where 

Q(a. 6, c) := cos (f (a — &)) [J(c, a, 6) + J(c, 6, a)] 
^ 3 ' 9 ^ - cos (f (a + 6)) [/(a, b, c) + J(a, b, c)] 
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and I(a, b, c), J(a, b, c) are the integrals defined by 

(3.10) I(a,b,c):= f C(l - a, q)((l - b, q)C(l - c, q) dq 

Jo 

and 

(3.11) J(a,b,c):= [ C(l-a,q)C(l-b,q)ai-c,l-q)dq, 

Jo 

where Q{z 1 q) is the Hurwitz zeta function. 

We consider the behavior of each of the factors in (|3. 8[) as b = e — ► 0. 

1) The term A(e) is regular in view of A(0) = 4^. The first few terms of its expansion 
are 

(3.12) A( £ ) = -L + (X±**M e + 0(e >). 

Z7T Z7T 

2) The expansion 

(3.13) cos (-(a ± e)J = cos — t - sin — e + 0(e 2 ) 
is elementary. 

3) To obtain the limiting value of Q(a,e,c) as e — > 0, it is required to analyze the 
behavior of ((s,q) as it appears in 

(3.14) I(a,e,c)= f C(o, <z)C(e, g)C(c, 9 ) dq. 

Jo 

The function C(e, q) has a pole at e = and its Laurent expansion at that pole is 

(3.15) C(e,q) = -- £ -m + 0(e). 

where ^{q) = logr(q). This expansion is not uniform in q, as the difference 

(3.16) C( £)(? ) -(-!/£-#?)) 

blows up as q — > 0, for any fixed e > 0. Indeed, 

(3.17) V(9) = -J-7 + 0(<z), 
and 

(3.18) ((e,q)=q e - 1 +t(l-e) + 0(q), 

showing that (|3.16p is not bounded as q — > 0. This issue is resolved by shifting the 
second argument of the integral in (|3.14p . Lemma |3~31 gives the integral J(a, e, c) in 
terms of integrals of the form I (a 1 , 1 + e, c'). The integrand here contains the factor 
C(l + e,q) = C(—e,q), whose expansion involves logL(<7), producing a similar blow 
up. Lemma 13.41 shifts the second argument again, and now I(a,e,c) is expressed in 
terms of I(a' , 2 + e, c'). The integrand now contains £(2 + e, q) = C(— 1 — e, 3) and 
a uniform expansion is finally achieved. 
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Lemma 3.3. For e > and a, c > 1 we have 

(3.19) I(a, e, c) = — [(a - l)J(a - 1, e + 1, c) + (c - l)J(a, £ + 1, c - 1)1 . 

e 

Similarly, 

(3.20) J{a,e,c) = -- [(a - l)J(a- l,e + l,c) - (c- 1) J(a, e + 1, c - 1)] . 
Proof. The identity 

(3.21) 

yields, for e > 0, 



C(e,g) = --f-C(l+e,«) 
e ag 



Z(a,e,c) 



1 



C(a, g)C(c, g)^-C(i + g) dg 

aq 



= -C(o,g)C(c,g)C(i + e,9) 



1 1 

o e 



C(l + e,g)— [C(o,g)C(c,g)] dg. 



Now observe that the boundary terms vanish for e > and a, c > 1, according to 
the identity (|3.5p . The identity for J follows along the same lines. □ 



The expansion of I (a, e, c) requires to iterate the result of Lemma 13.31 twice. 
Lemma 3.4. For e > and a, c > 2, 

I{a, e, c) = ^jyy [ (a ~ ~ 2)/(a ~ 2 ' £ + 2 ' c) 

+ 2(a - l)(c - l)J(a - 1, e + 2, c - 1) + (c - l)(c - 2)J(o, £ + 2, c - 2)] , 
with similar expressions for J(a, e, c) and J(c, a, e). 
Replacing in the expression (|3.8p gives 

(3.22) T(a, £, c) = «^ sin ( - ) LT(a, e, c) 
where 

(3.23) LT(a, e, c) = cos (^( a ~ e )) ^+( a : e i c ) — cos (^( a + e )) ^-( a ' e > c )i 

with 
(3.24) 

iJ+(a,e,c) = (1 - a) (2 - a) J(c, a - 2, e + 2) + 2(1 - a)(l - c)J(c- l,o - 1,£ + 2) 
+ (1 - c)(2 - c) J(c - 2, o, e + 2) + (1 - c)(2 - c) J(c - 2, e + 2, a) 

- 2(1 - a)(l - c) J(c - 1, e + 2, a - 1) + (1 - a) (2 - a) J(c, e + 2, a - 2), 

and 
(3.25) 

H-(a,e,c) = (l-o)(2- a)J(o- 2, e + 2, c) +2(1- a)(l - c)J(a- l,e + 2,c- 1) 
+ (1 - c)(2 - c)J(a, £ + 2, c - 2) + (1 - a)(2 - a) J(a - 2, e + 2, c) 

- 2(1 - a)(l - c) J(a - 1, £ + 2, c - 1) + (1 - c)(2 - c) J(a, e + 2, c - 2). 
The last factor in (|3.22[) is now expanded in powers of e to obtain 

(3.26) LT(a,e, c) = C a (M+-M_)+e [c a (N + - N_) + |s a (M + + M_)] +0(e 2 ), 
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with a self-explanatory notation. For example M+ is the limit as e — > of the 
expression multiplying cos(§(a — e)) in (|3.23|) . N+ is the term of order e in the 
same factor and C a = cos(-7ra/2). We now reduce all the terms in (|3.26p . 

Proposition 3.5. The identity M + = M_ holds. Therefore, the term LT(a,e) is 
of order e, and we conclude that T(a, e, c), in (|3.22[) , is non-singular as e — > 0. 

Proof. The proof of this result begins with 

(3.27) 

M+ = (1 - a) (2 - o)J(c, a - 2, 2) + 2(1 - a)(l - c) J(c - 1, a -1,2) 

+ (1 - c)(2 - c) J(c - 2, a, 2) + +(1 - c)(2 - c) J(c - 2, 2, a) 

- 2(1 - a)(l - c) J(c - 1, 2, a - 1) + (1 - a) (2 - a) J(c, 2, a - 2), 

and 

M_ = (1 - a) (2 - a)J(a - 2, 2, c) + 2(1 - a)(l - c)J(a - 1, 2, c-1) 

+ (1 - c)(2 - c)J(a, 2, c - 2) + (1 - a) (2 - a) J(a - 2, 2, c) 

- 2(1 - a)(l - c) J (a - 1, 2, c - 1) + (1 - c)(2 - c) J(a, 2, c - 2). 

We now establish some relations for the integrals I and J. These are then 
employed to show that M + = M_ . 

Lemma 3.6. The integrals I and J satisfy 

(3.28) J(z, z', 2) = I{z, z', 2) and J(z, 2, z') = J(z', 2, z). 
Proof. Start with 

(3.29) J(2,2',2)= / C(«.«)C(*'.9)C(2,i-?)dg, 

Jo 

with ((2, q) — ((— 1, g). From C(l — k,q) = —Bk(q)/k and the symmetry of B2 we 
obtain 

(3.30) C(2, 1 - q) = C(~l, 1 - 9) - C(~l, «) = C(2, g). 

The first identity follows from there. A similar argument establishes the second 
one. □ 

A direct calculation using Lemma \5M now shows that M+ — M_ in (|3.26|) . □ 
Calculation of M+. The next step is to simplify the expression for M + given in 

(ED . 

Proposition 3.7. Let a > 2, c> 2 e 1. Then 

C(a + c)B(a, c) 



(3.31) M+ = 2 C(«)C(c 



1 — tan(7ra/2) tan(7rc/2) 
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Proof. Start by producing an alternative form of the term J(c, a — 2, 2) given by 

J(c,a- 2,2) = [ C(c,q)C(a-2,q)a2,l-q)dq 
Jo 

J C(c,g)C(2,l-g)^C(o-l,g)do 
C(a-l,q)^- [C(c,q)C(2,l-q)] dq 
C(o-l,g) [(c-l)C(c-l,«)C(2,l- g ) 

+C(c,?)C(l,l-?)] 

C(i,«) = -C(i,i-«), 



1 


ft 




2 




1 




2 




a 




1 




2 




a 



Now recall that 
since f (1, q) = —Bi(q) = | — g, and also 



C(2,5)=C(2,l-«), 
as shown before. Conclude that 

(3.32) (a - 2) J(c, a - 2) = (1 - c)I{a - 1, c - 1, 2) + /(a - 1, c, 1). 

Reversing the order of a and c and considering the expression in (|3.27p . the first 
three terms of M + in (|3.27p are reduced to 

(3.33) (1 - a)I{a - 1, c, 1) + (1 - c)I(c - 1, a, 1). 

A similar analysis gives an expression for the last three terms in (|3.27p . This yields 
a formula for M+ in terms of the integrals / and J where one of the parameters is 
1. 

Lemma 3.8. The term M+ in (|3.27|) is given by 

A/+ = (1-0)7(0-1, c, l)+(l-c)/(c-l, a, l)+(l-o) J(o-l, 1, c)-(l-c) J(a, l,c-l). 

We now consider the integrals appearing in this expression for M+. 
Lemma 3.9. The integral J satisfies 

(3.34) J(u,l,v) = -J(v,l,u), 

Proof. This comes directly from £(1, 1 — q) = — £ (1, q). □ 

The expression for M+ in Lemma 13.81 is now simplified. 
Lemma 3.10. Let a > 2, c> 2 e M. Then 

(1 - o)7(o - 1, c, 1) + (1 - c)7(o, c - 1, 1) = C(«)C(c) - / C(o. ?)C(c, 9) ^ 



and 

(1 - a)J{a - 1, 1, c) - (1 - c) J(o, 1, c - 1) = C(o)C(c) - f C(a, q)((c, 1 - g) dry. 

Jo 

Therefore 

(3.35) M+ = 2C(o)C(c) - / C(a,q)((c,q)dq- [ C( a , g)C(c, 1 - g) dg. 
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Proof. Start with 

(1 - a)I(a - l,c, 1) = (1 - a) I ((a - 1, <?)C(c, q)((l, q) dq 



and the identity (a — l)C( a — 1, Q 1 ) = ^C( a i Integration by parts and the value 
£(1, q) — \ — q give the first identity. The rest of the formulas are established in a 
similar form. □ 

The identites 

,cos(|(a-c)) 



f C(l - a, q)C{l - c, q) dq = ^{l-a- c)B(a, 
Jo 



c 



cos(f (a + c)) 



and 

fi 



C(l - a, g)C(l - c, 1 - q) dq = C(l - a - c)B(o, c) 

appear in Theorem 3.1 in [4]. Replacing in (|3.35[) produces the final expression for 
M+ claimed in (|3~3T|) . □ 



Calculation of N + . This is defined in (|3.26[) as the term of order e in the expansion 
of H + (a, e, c) defined in (|3.24[) . We illustrate the general method of calculation by 
computing the term of order e in the integral: 

J(c, a - 2, e + 2) = f ((c, q)((a - 2, q)((-l -e,l-q) dq. 
Jo 

The term of order s is 

C(c,q)C(a-2,q)C(-l,l-q)dq 



1 



a- 2 



— C(o-l,g) )((c,q)C(-l,l-q) dq. 



Integrate by parts and observe that the boundary terms vanish to see that the term 
of order e is 

f aa-l,q)^-[ac,q)a-l,l-q)] dq. 
a - 2 J dq 

(a — 2) x the term of order e in J (a — 2, c, e + 2) 



Thus, 



is 



(c - 1) / C(c - 1, q)((a - 1, <?)<'(-!, 1 - (?) dq 



C(c,g)C(o-l,«)^C'(-l,«) 



u—l — q 



dq. 



Now use 

| I CVM) = ^L = _ 1 | I C(^ W ) = ^L = _ l (-<(, + M)) = -C(0,.) + C'(0 !W ) ) 
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to check that the contribution coming from 1, 1 — q) vanishes and to verify the 
result described next. 

Lemma 3.11. The part of N+ coming from integrals with e+2 in the third variable 
is 



C(c, q)((a - 1, q)C(0, 1 - q) dq - / ((c, q)((a - 1, q)£(0, 1 - q) dq 



C(o, q)((c - 1, q)C(0, l-q)dq- / £(<»> q)((c - 1, g)C'(0, 1 - «) dq 



-(1-a) 



-(1-c) 



Similar calculations produce the other parts of 7V + . We spare the reader the 
details. 

Proposition 3.12. The term N + is given by 

2C(o)C(c)-/ ((a,q)((c,q)dq- f ((a,l - q)((c,q) dq 



+ (l-o) 



+ (1"C) 



C(c,g)C(a-l,5)C'(0,l-<z) dg 



+ / C(c,i- 9 )C(a-i,g)C'(0,i-g)dg 



C(a,g)C(c-l,g)C'(0,l-g) dg 



C(a,l- ? )C(c-l,«)C / (O s? )dg 



The process is now repeated to produce a similar expression for 7V_ . The result 
is stated next. 

Proposition 3.13. The difference N + — N- is given by 

N+ - N- = (1 - a) [ C(a - 1, ?)C(c, g) [C'(0, g) + C'(0, 1 - q)} dq 



+ (1 - a) / C(a - 1, 9)C(c, 1 - g) [C(0, g) + C'(0, 1 - g)] dg 
Jo 

+ (l - c) / c(o, g)C(c - l, g) [C'(0, g) + C(o, l - g)] dg 



(l - c) / C(o, g)C(c - 1, 1 - g) [C'(0, g) + C'(0, l - g)] dq. 



The difference N + — iV_ is now given in the form stated in Theorem [2TTJ Recall 
that that ('(0, g) = logT(g) - logv 7 ^, so that 

(3.36) C'(0, q) + C'(0, 1 - q) = log [T{q)Y{l - q)) - log27r = - log(2 sin^q), 
where we have used the reflection rule for the Gamma function 

7T 



(3.37) 

Integration by parts gives 



T{q)T{l-q) = 



sm irq 



(1-a) ((a - 1, q)C(c, q) dq = -(1 - c) / C(a> «X(c - 1, g) dg, 
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and 

(1 - a) [ ((a - 1, q)C(c, 1 - q) dq = (1 - c) / ((a, q)((c -1,1-q) dq. 
Jo Jo 

These identities show that in the expansion of N + — N- in Proposition 13 . 1 31 we 

can ignore the contribution of the q-independent term — log 2 in (|3.36p and produce 

our final expression for the term N+ — iV_ . 

Proposition 3.14. The term N+ — N- in (|3.26|) is given by 

^ [C(c, q) (({a, q) - ((a, 1 - q))] logsin(Trq) dq 

l 



= -7T / C(c, q) [((a, q) - ({a, 1 - q)] cot{nq) dq. 
Jo 

Taking the limit as e — > in (|3.22p , we obtain the result described in Theorem 



4. An expression for the Tornheim sum T(m, 0, n) 

In this section we analyze the behavior of T(o, 0, c) given in Theorem 12. II as the 
parameters a and c approach positive integer values. The notation a = m + e a , c = 
n + e c with m, n € N, m, n > 2 and e a , e a — ► 0, is empoyed throughout. 

Observe that X(z) is singular as z becomes a positive integer, so in order to 
obtain the limiting value T(m, 0, n), it remains to consider the limiting behavior of 
T(a,0,c). 

Define 



(4.1) Li(m, n) := lim lim 4A(a)A(c) sin ( — ) sin 

a — >m c — >n V 2 / 



x U (a)C (c) - 



1 - tan (^f ) tan (^) 



and 

(4.2) L 2 (m, n) := - lim lim 2A(a)A(e) cos ( — sin ( — 

x / [C(a, q) ~ ((a, 1 -?)] C(c,q) cot 7rgd<j. 



The expressions for Li(m,n) and L^ivn^n) given in this section give Theorem 12.21 
from Theorem 12. II 



The evaluation of L\(m,n). This is elementary. The simplification employs the 
relations 

(4.3) ^ 2fc ) = ( " 1) Tff fcg2fc ' keNU{0}, 

/ 1 \fc+l D 

C(l-fc)= l j fc , fceN, 
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and 

Proposition 4.1. The function L 1 (m, n) defined in (|4.f [) is given by 

(4.4) Li(m 7 n) = £(m)£(n) — — C( m + 
Proof. The expansion of A(z) about z = A; G N is 

(4.5) A(* + e) = (-l)""" 



i + log 2?r - -0 (fc) + O (e) 



T(fc) 

We now examine the behavior of the other factors in (14. ip as e a ,e c tend to zero. 
The result depends on the parities of to and n. 

To simplify our notation we define the A-parity of the integer fc, pk, as 

(-l) fe / 2 , k even, 



(4.6) Pfc :=(-l)Lk/2J = 

Then, for k even: 



(_l)(*-l)/2 j jfc odd) 



(4.7) sm(^(k + e))=^p k e + 0(e 2 ), 

(4.8) tan(|(fc + e)) =| £ + 0(e 2 ), 

(4.9) C(fc + e) = C(l-fc) + 0(e). 
and for fc odd: 

(4.10) sin(|(fc + e)) =p fc + 0(e), 

(4.11) tan(J(fc + e)) = -— + O(e ), 

\ 2 / 7T£ 



(4.12) C(fc + e) = -C'(l-fc)e + 0(e 2 ), fc > 3, 

since C(-2fc) = for fc E N. 

Therefore 

(a) m and n even: 

(27r) m+ ™ 

Li(m,n) = . [C (1 - m) C (1 - n) - B(m,n)C (1 - m - n)]p TO p„ 

41 (to) 1 (nj 

(b) to even and n odd: 

f27r) m ^™ — 1 

ii(m,n) = [C(l - ™) C (1 -n) - B(m,n)C (1 -m-n)]p m p„ 

1 (to) 1 (n) 
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(c) m odd and n even: 
(27r) m ^™~ 1 

L 1 (m,n) = [C (1 - to) C (1 - n) - B(m, n)(' (1 - m - n)] p m p r , 

1 (to) 1 (n) 



(d) to and n odd: 

\m-\-n 

4r (to) r ( n ) 

The result now follows by using the relations (|4.3[) . □ 



T ( \ ^ 

L\ (to, n) 



4 

-~ C (1 - to) C (1 - n) + B(m, n)C (1 - to - n) 



n ■ 



The evaluation of L2(m,n). This proceeds along similar lines. We employ the 
expansion 

C(k + e,q) = ((k,q)+e('(k,q) + o(e) 

= -^[B k {q)+eA k (q)]+o{e), 

and the reflection property of the Bernoulli polynomials, 
(4.13) B k (l-q) = (-l) k B k (q), 

to establish the vanishing of some of the integrals that appear in intermediate calcu- 
lations. The results are expressed in terms of the four integrals Iaa, Iab, Ibb, Jaa 
introduced in (12.51). 



Proposition 4.2. The limit L2(m,n) defined in (|4.2p is given by 

(27r) m+ra_1 

(4.14) L 2 (m,n)=p m+n — £ 2 (m,n). 

mini 

Here p n + m is defined in (|4.6p and the function li{m, n) is given in terms of the 
following basic integrals, 

(4.15) I BB (k,l):= f B k (q)B l (q)K(q)dq, 

Jo 

i r 1 

(4.16) I AB {k,l):=- A k {q)B l {q)K{q)dq, 



1 ^ 



(4.17) lAA{k,l) := s A k (q)Mq)K(q)dq, 



1 f 1 

(4.18) J AA {k, A k {q)A[(l - q)K(q) dq. 

n Jo 

by 

(a) to and n even: 

(4.19) l-x(jfit n) = mlABi'm — 1,71.) + i/^fra, n — 1), 

(b) to and n odd: 

(4.20) n) = mlAB{n, m — 1) + nlABin — 1, to). 
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(c) m odd and n even: 

(4.21) £ 2 (m,n) = - (mI B B(m — l,n) + nI B B{m, n—1)), 

(d) m even and n odd: 

(4.22) £ 2 (W) = — mlAA^jn — 1) 7i) — nlAA(m, n—1) 

— mJAA{m — l,n) + nJ AA(jn, n — 1), 
Proof. We make use of the expansion 

/7T \ fpfe + o(e), A; even 

(4.23) cos( J(fc + e)) = T tt ' . , 

V ; V2 V V j^-pfe - e + ( e ), fe odd. 

In the case to even, the singularity of A(m + e a ) as e a — ^ is balanced by 

C (m + E a , q) - C (m + £„, 1 - g) = -— L4 m (<?) - A m (l - q)] e a + O (e a 2 ) . 

m 

For m odd, the combination A(m + e ) cos ((m + e a )7r/2) is regular as e a 0. 
The singularity at n + s c is treated similarly: the term A(n + e c ) sin ((n + e c )Tr/2) 
is regular for n even; for n odd, the singularity of X(n + e c ) requires the vanishing 
of the integral in (|4.2p . This is a consequence of the symmetry of the function 
[C(a,q) — ((a, 1 — q)] cot7rg, about q — 1/2. 

Introduce the notation 

(27r) m+n - 2 



(4.24) C^m.n ■ — Pm-\-n~ 

where p m +n is defined in (|4.6[) . Then 



>!nl 



(a) for to and n even: 

L 2 (m,n) = -na m , n / [A m (g) - A m (l - q)]B n (q) cotwqdq, 
Jo 

(b) for to even and n odd: 

L 2 (to, n) = 2a mj „ / [A m (q) - A m (l - q)]A n (q) cot wqdq, 



(c) for to odd and n even: 

L 2 (m, ra) = -n 2 a m . n / B m (q)B n (q) cot irqdq, 
Jo 

(d) for to and n odd: 

L 2 (m,n) = -27ra TO ,„ / S m (g)A n (g) cotTrgdg. 



The final expression in Proposition 14.21 is obtained by writing the integrals above 
in terms of the basic integrals defined in (|4. 15|) . This can be achieved by using 

tt cot nq = — log sin nq = K'(q), 
dq 
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integrating by parts and employing the recurrence relations 

(4.25) —B k {q) = kB k ^{q), 

(4.26) 4-Mo) = kAfc-ifa) + --^-Bk^q), k>2. 
dq K — l 

See [5] for a derivation of (|4.26[) . 

The kernel K(q) has only a logarithmic singularity at q = 0, so it is possible to 
separate the integrals involving the difference A k {q) — A k (l — q). The functions 
Ak(q) are well behaved in [0, 1]. Only A 1 (q) = log [r(g)/v27n has a singularity at 
q = 0, and this is only logarithmic, and therefore integrable in the cases of interest 
here. □ 

Example 4.1. The Tornheim sum T(2, 0, 2) is the simplest example of even weight. 
Its exact value is 

(4-27) T (2,0,2) = |c(4) = ^, 

see Theorem [231 states that T(2, 0, 2) = Li(2, 2) + L 2 (2, 2), with 

(4.28) L!(2,2) = ^, 
and 

(4.29) L 2 (2,2) =47r 3 (/ AB (l,2)+/ AS (2,l)). 
Therefore 

(4.30) T(2, 0, 2) = ^ + 4tt 2 jf* A 1 {q)B 2 {q)K{q) dq + in 2 J* A 2 (q)B 1 (q)K(q) dq, 
and we have the evaluation 

(4.31) jf [A 1 (q)B 2 (q) + A 2 (q)B 1 (q)} log sin irqdq = — . 
Example 4.2. Theorem 12.21 gives 



(4.32) T(2, 0, 6) = ^ + ^/ab(1, 6) + ^ab(2, 5). 

Theorem 4.3. Every Tornheim sum can be expressed as a finite sum of integrals 
of the form Iab- 



5. A REPRESENTATION FOR /,4b (m,n) IN TERMS OF CLAUSEN FUNCTIONS 

In this section we consider the integral 



1 



(5.1) I A B{m,n) = - A m (q)B n (q)K{q) dq, 

n Jo 

and express it in terms of the family introduced in (|2.11|) . 

(5.2) X Kl ■= (-l)L'/aJ JL [ log T(q)B k (q) Cl l+1 (2irq) dq. 

Jo 

Here CI/ (x) are the Clausen functions defined in (|1.18p and (|1 . . 
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The reduction of the integrals Iab employs a new family of special functions, 
K n (q), n £ No, defined as the iterated primitives of the kernel 

(5.3) K (q) := -K(q) = -logsinTrq, 
that is, 

(5.4) K n (q) n f K n ^{q')dq', n > 1. 

Jo 

The functions K n (q) are normalized by the condition K n (Q) = for n > 1. 

Note 5.1. The plan is to use (n + l)K n (q) = K' n+1 (q) and integrate by parts to 
transform the expression for lAg(m, n) into a finite sum of integrals in which the 
only type-A function that appears is Ai(q) ~ logr(q) — log-v^r. 

The Fourier expansion of Ko(q) is 

(5.5) K (q) = lo^ + icir (27T9), 

where Cli is a Clausen function. We first show that K n (q) is the sum of a polynomial 
in q of degree n and a multiple of the Clausen function C\ n +i{2irq) (see (|5.22[) ). We 
then show that the polynomial part can be integrated out explicitly in terms of zcta 
values, leading finally to the analytic expression for the Tornheim sums T(m, 0,n) 
of even weight N — m + n in terms of the integrals Xki, given in Thcorcm l2.3l 



We first recall some basic properties of the functions A m (q), introduced in [5]. 
The function 

d 

(5.6) 4(9 : = m TrC(z,q) 

az 

satisfies, for to G N — {1}, the recurrence 

(5.7) -rA m {q) = m4-i(«) + -B m _ 1 (g), 

aq to — 1 

with initial condition 

(5.8) A 1 {q) =logr(< 7 )-lo g v / 2^. 
The boundary values are given by 

(5.9) A m (0) = A m (l) = toC'(1 - to), for to > 2, 
and 

(5.10) A x (0+) = -oo, Ai (1) = - log v 7 ^. 
The recurrence (15.71) is similar to the relation 



(5.11) -^-B m (q) = mB m -i(q), 

aq 

satisfied by the Bernoulli polynomials. 

Theorem 5.2. The integral lAB(rn,n) can be expressed in terms of the family 

(5.12) Uk t i := [ \o g r(q)B k (q)K l (q)dq, 

Jo 

with k + I = to + n — 1. 
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Proof. The proof is by induction on m. We show first that any integral of the form 



f A m (q)B n (q)K p (q)dq 
Jo 



can be expressed in terms of 
(5.13) U*, := / A l {q)B k {q)K l {q)dq, 



with k + I + 1 = in + n + p, and then transform into Uk,i by using A±(q) 
lo g r(g)-logV2^. 

Start with 

(p+1) / A m {q)B n (q)K p (q) dq = [ A m (q)B n (q)-^-K p+1 (q) dq. 

Jo Jo d q 

Integrate by parts and use K p+ i(0) = 0, (|5.7|) . (|5.9|) and (|5.1ip to derive 



(5.14) 



(p+1) / A ro ( 9 ) J B n (g)iL p ( 9 ) = A m (l)B n (l)K p+1 (l) 

m I A m _ 1 (q)B n (q)K p+1 (q) dq 

B m ^ 1 (q)B n (q)K p+1 (q) dq 



o 



m — 1 

-n/ J 4 m (g)S n _ 1 (g)isT p+ i(g). 
Jo 

The term involving the product of two Bernoulli polynomials can be evaluated 
in closed form, and only the last term in (]5 . 14[) requires further analysis. Repeating 
the process to this last term produces 

(p+l)(p + 2) / A m (q)B n - 1 (q)K p+1 {q)dq = A m (l)B n - 1 {l)K p+2 (l) 
Jo 

-ml A m -i(q)B n -i(q)Kp +2 (q)dq -/ B m _ 1 (q)B n _ 1 (q)K p+2 (q) dq 

Jo m ^ 1 Jo 

-(n-1) A m (q)B n - 2 {q)K p+2 (q)dq. 
Jo 

We conclude that each integration by parts produces terms that can be evaluated 
in the stated closed form (by the induction hypothesis) and one integral where the 
index of the Bernoulli polynomial term is decreased by one. Note that the indices 
j, k, I of the three functions in any term of the integrand appearing in this procedure 
satisfy j + k + l = m + n+p. The case of interest, /^(m, n), corresponds to p = 0. 
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Eventually we arrive at 
(p + n + 1) / A m (q)B (q)K p+n (q)dq = A m (l)B (l)K p+n+1 (l) 



' ^p+n+i (<z)-Bo (g) I mA m -\ (q) H — --B m _i(g) ] 

o V m- i y 



A n (9)^--B (g)i^ p+ „ +1 (g) dq. 

/o 

The fact that B (q) = 1, so that the last integral vanishes identically, completes 
the argument. 

It remains to prove that integrals of the form 



(5.15) 



V, 



B k (q)B m {q)K n (q)dq 



can be evaluated in closed form and to use the relation (|5.8|) to eliminate Ai(q) 
from the formulas. 

The first step in this part of the proof is to use the relation 



fe(ni,n 2 ) r 

(5.16) B ni (q)B n M= E 

k=0 



n-2 



n ^2k +n \2k 



Til 



B 2 k 



n\ + n 2 — 2k 



B ni+n2 -2k(q) 



(ni + n 2 )\ 



with 
(5.17) 

for the product of two Bernoulli polynomials in order to reduce Vk, m ,n to sums of 
terms of the form 



fc(m,n2) =Max{ln 1 /2j}, [n 2 /2\}} , 



(5.18) 



W, 



, B m (q)K n (q)dq. 
o 



Integration by parts show that (|5.18[) can be written in terms of the M^-integrals 
with second index 0. Indeed, 



W m ,n = 



If 1 d 

—r / K n (q)—B m+1 (q)dq 
+ 1 Jo "9 



B m+1 {\)K n {\) 1 



to + 1 to + 

B m+1 {\)K n {l) n 



771+1 



m + 1 

and iterating this procedure yields 

(5.19) W m ,n = V t — J -—B m+k (l)K n - k+1 (l) + 

f— ' TO + 1 )k 



f 1 d 
- I B m+1 (q)—K n {q) 
1 Jo dq 

W^m+l.n-l, 



fe=l 



(-1)" 
(TO + 1), 



where (m)k — m(m +!)■■■ (m + fc — 1) is the Pochhammer symbol. 
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The closed form of Vk, m ,n follows from 

i (log 2 r = 0, 

(5.20) W rfi = j B r (q)K (q)dq= (0 r odd, 

[(-l) r / 2 r!C(r + l)/(27r) r , r > even, 

given as Example 5.2 in 0]. □ 

Note 5.3. The boundary terms K n (l) are given in (|5 . 29[) . 

Our final step in the reduction of the integrals Tab makes use of a representation 
of the kernels K n (q) in terms of the Clausen functions Cl„(x), defined in (|1.18[) and 
Oi . 

Proposition 5.4. The kernels K n (q) are given by 

(5.21) K (q) = ]og2 + ±CL 1 {2*q), 
and, for n > 1, 

(5.22) ^)= g "log2 + n! £ \J )2fc( ^_ 2 fc)! ^ + P " (2^ ^M, 
where the coefficient p n = (— 1)L"/ 2 J j s the A-parity of n defined in (|4.6|) . 

Proof. The Fourier series expansion of the function Ko(q) — — logsin(7rg) 

«„s » . , „ 1 \ -> cos(27rfco) 

(5.23) K (q) =log2 + -^ V fc ^ . 

fc=i 

is standard. To compute the expansion for Ki(q), we cannot just integrate term by 
term the series (15.23[) . since it is not uniformly convergent for q £ [0,1]. To bypass 
this problem, observe that the indefinite integral of Ko(q) can be written as 

(5.24) J K ( q )dq = q\og (l - e 2 ^) - q log sin (irq) - ^q 2 - ± Li 2 {<£*««) , 
where Li n (z) is the polylogarithmic function defined for z E C, \z\ < 1 by 

OO 

(5.25) U n(z)'-=J2]^> n ^ 2 - 

k—l 

Separating the real and imaginary parts we obtain, for < q < 1, 

sin(27rfc(7) 



k 2 

k=l 



f 1 

(5.26) K (q)dq = qlog2+ — Y^ 

^ 71 k=l 

7T . ,1 cos(27rA:(7) 
fe=l 

and the proposition is proved for n = 1 since the term in brackets is a ^-independent 
constant. This can be seen directly from Hurwitz's Fourier series representation of 
the Bernoulli polynomials 
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where the prime indicates that the term k = must be excluded in the sum. 
The result for n > 2 can be obtained directly from the expression for K±(q), 

,~ „„s j . , 1 sin(27rfco) 

(5.28) KM = g log2 + - \ 2 * 

fc=i 

integrating successively term by term, which is now legitimate since, for n > 2, the 
series defining the Clausen function Cl„(x) is uniformly convergent. □ 

Note 5.5. The evaluation of the basic integral Iab{iti, n) requires the boundary 
values of K n (q), n G N, at q = 1. These are given by 

L(«-i)/2j f _ lV fc+Vf2A- + -n 

(5.29) «„(!> = l„ g 2 + „! g y^gg^ . «N. 

where we have used the special values Ci2n(27r) = and Cl2 n +i(27r) = ((2n +1). 

Note 5.6. The method presented in this section reduces the evaluation of all even- 
weight Tornheim sums to the evaluation of the integrals U m . n , defined in (|5.12[) . 
The expansion (|5.22|) and the formula 

n+l . 

,1 1 L 2 J 



(5.30) I fWM*,-^ £ (-l) fe (; + _ 1 1 )J§2T[^)-C'(2^] 

- — ^ 2fc J^p« 2fc+1 ) + 



n+l 



k=l 

(with 5 = 2 log \/2tt+j, where 7 is Euler constant) given as (6.14) in [4], reduces the 
evaluation of Tornheim sums to the evaluation of the integrals X mj „, introduced in 
(|5.2p . The result described in Theorem 12.31 appears from transforming the explicit 
representation (|2.6|) in Theorem [221 by the methods presented in this section. The 
details are left to the reader. 

Note 5.7. The value 

(5.31) x ,„ = (-l)L»/2J^L f\ogT( q )C\ n+1 (2TTq)dq 

[ Z7r > Jo 

can be obtained from the integrals 

A + log k 



[ log T(q) sm(2Trkq) dq = 

f 1 1 

J log T(q) cos(27rfcg) dq = — 



where A = log(27r) + 7, with 7 the Euler constant. This appears in [5] 6.443.1 and 
6.443.3. It follows that 



(5.32) X ,„ = (-1)L"/ 2 J 



nl 



J J^ F (AC(n + 2)-C / (n + l)), for n odd, 
(2tt)" + for 7i even. 



For our evaluations of even- weight Tornheim sums, we will only need the result for 
n even. 
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Example 5.8. This is a continuation of Example 14.11 We express here the sum 
T(2, 0, 2) first in terms of the [/-integrals. Theorem 12. 2\ with m = n = 2, gives 



T(2,0,2) = C(2) 2 -ic(4)-47r 2 



A 1 (q)B 2 (q)K a (q)dq 



A 2 (q)B 1 (q)K a {q)dq 



o 



The second integral is now modified to avoid the presence of the function A 2 (q) . 
Write K (q) = (d/dq)Ki(q), integrate by parts and use formulas (|5.7[) and (|5.11|) 
to get 



A 2 (q)B 1 (q)K (q)dq = A 2 (l)B 1 (l)K 1 (l)- / A 2 (q)B (q)K 1 (q)dq 

Jo 

-2 f A 1 (q)B 1 (q)K 1 (q)dq-2 f Bl{q)K x {q)dq. 
Jo Jo 

The new integral involving A 2 (q)Bo(q)Ki(q) can be dealt with in a similar manner, 
to produce 

jf A 2 (q)B 1 (q)K (q) dq = A 2 {l)B 1 {l)K 1 {\) - -A 2 {\)B a (\)K 2 {\) 

+ f A 1 (q)B a (q)K 2 (q)dq-2 [ A l {q)B 1 {q)K 1 (q) dq 
Jo Jo 

B (q)B 1 (q)K 2 (q)dq~2 / Bf(q)Ki(q) dq. 
) Jo 



From (|5.9p and (|5 . 29[) we see that the boundary terms calcel out and therefore the 
Tornhcim sum T(2, 0, 2) is given by 

T(2, 0, 2) = C 2 (2) - C(4)/2 - 4vr 2 (t/* - 2^ + U* <2 + Vfc.i.a - 2V W ) . 

The expression for Vk, m ,n shows that Vo,i,2 = 214,1, l = jo log 2. Therefore, the sum 
T(2, 0, 2) is expressed in terms of the ?7*-integrals as 

T{2, 0, 2) = ^ - 4tt 2 (t/* Q - 2C/* X + £/ * 2 ) . 
The transformation from n to £/ mj „ gives 

T(2, 0, 2) = ^ + 1 log(27r) (tt 2 log(2) + 9C(3)) - 4tt 2 (f/ 2>0 - 2[/ M + L/b )2 ) . 
The final transformation is to write U m , n in terms of X m ^ n to obtain 

(5-33) T(2,0,2) = ^-y 2 * 2 , 

where 

(5.34) y 2 * 2 := 4tt 2 (X 0>2 - 2X 1A + X 2fi ) - 2((3) log(27r). 

Example 5.9. At weight 8, we consider first the value of 

7T 8 47T 7 

T(2, 0, 6) = — + — (2I AB (1, 6) + 6I AB (2, 5)) . 
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Expressing the Iab integrals in terms of the U* integrals we find 
2J ab (1,6) + 6/ab(2,5) = 

- - (US fi - 6C/r.5 + 15^,4 - 20£/ 3 * 3 + 15C/ 4 * 2 - GUI, + Ko) , 

7T 

and in terms of the U integrals: 

9T nn + KT OK 1 W9^f log(2) C(3) !5C(5) 315C(7) \ 
2I AB {1, 6) + QI AB {2, 5) = - log( 27 r) ^— - — - + -^-j- j 

- - (U as - 6f/i, 5 + 15[/ 2 , 4 - 20C/ 3 . 3 + 15(7 4 . 2 - 6(7 5 .i + tf 6 ,o) • 

7T 

Thus we obtain the representation 

, r „„. m/ „ „ ^ 7T 8 27T 6 , ,„ . /log(2) C(3) 15C(5) 315C(7)\ 

5.35) T (2,0,6) = h log(2?r) 6W - 2±-L 2±-L -\ 

v ' v ' ' ; 8100 45 sv ; V 21 2?r 27r 27r / 

87T 6 

- -7=- (U fi - 6U ll6 + 15U 2A - 20£/ 3 , 3 + 15U 4 ,2 - 6U 5 ,i + U 6 , ) . 
45 

In terms of the X m>n integrals we have 

(5-36) T (2,0,6) = -^-Y 6 * 2 , 

where 

(5.37) F 6 * 2 := ^- (X 0fi - 6X1,5 + 15X 2j4 - 20X 3 , 3 + 15X 4 , 2 - 6X5,1 + X 6 , ) 

-6C(7) Iog(27r). 

Example 5.10. The other two Tornhcim sums of weight 8 are T(3,0, 5) and 
T(4,0,4). They are given by 

r(3, 0, 5) = + C(3)C(5) + ^ (5J AB (4, 3) + 3I AB (5, 2)) 

and 

T(4 ' °' 4) = 1^5 + T~ {Iab{3 > 4) + /AS (4 ' 3)) ' 

respectively. 

Expressing the Iab integrals in terms of the X integrals we find 

(5.38) T(3, 0, 5) = - 6C(3)C(5) - F 3 * 5 , 
where 

S27T 6 

(5.39) r 3 * 5 := — (X ,6 - 3X1,5 + 3X 2 , 4 - X 3 , 3 ) - 15C(7) log(27r), 
and 

(5.40) T(4, 0, 4) = ^ - 4C(3)C(5) - Y 4 * A , 
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where 

87T 6 

(5.41) Y 4 * 4 := — (JC ,6 - 4Xi j8 + 6A 2 , 4 - 4X 3 , 3 + X 4 . 2 ) - 20C(7) log(27r). 

Note 5.11. In the examples given above, the integrals X^,i appear in a symmetric 
form. This is a general feature, as stated in Thcorcm l2.31 Depending on the parity 
of the integers m and n, the even-weight Tornheim sum T(m, 0, n) contains either 
Y m>n or Y n>m , where 

Ym - n : " m!(n-2)! ^' 1 : ' ( y j A ' " 2 /• 

There are linear relations among the integral Y mn of fixed weight N := m + n. 
These come from the linear relations among the Tornheim sums T(m, 0, n) of the 
same weight, discussed in the next section. For example, the identity 

T(2,0,6) + T(6,0,2) = ^C(8) 

gives 

Y 2fi + Y 6<2 = 12 log(27r)C(7) + |c(8) - 6C(3)C(5), 

5T(6, 0, 2) + 2T(5, 0, 3) = — C(8) - 8C(3)C(5) 

29 

5Y 2fi + 2r 5 , 3 = 601og(27r)C(7) + -((8) - 30C(3)C(5). 
A systematic study of the integrals Y m>n will be presented elsewhere. 



whereas 



gives 



6. A SYSTEMATIC LIST OF EXAMPLES 

Here we present a systematic evaluation of the Tornheim sums T(to, k, n) with 
to, k, n £ NU {0}. The sums are organized according to the weight N = m + k + n. 
The conditions to + n > 2, k + n > 2 and N > 3 are imposed for convergence. The 
symmetry relation T(m, k, n) = T(k, to, n) is used to impose m > k. 

We use the notation Zjq and Z® N introduced in the introduction. Huard's result 
(|1.4p allows us to evaluate any Tornheim sum in Z^ in terms of the sums in Z^. 
Before detailing a systematic algorithm to evaluate all the sums in Z^ for a given 
weight N and giving specific examples, we shall determine how many of the Torn- 
heim sums in Z^, for a given even weight N, remain undetermined after using all 
the known algebraic identities for these sums. 

As we discussed in the introduction, two of the sums in Z^ have known explicit 
evaluations: 



(6.1) 
and 

(6.2) 



T(0,0,N) = ((N-1)-((N), N>3, 



r(i,o,j\r-i) = 



N-2 



(N-1)C(N)- 53C(i)C(JV-0 



i=2 



TV > 3. 
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Euler proved that for m > 2, n > 2 the sums T(m, 0, n) satisfy the symmetrized 
identity 

(6.3) T(m, 0, n) + T{n, 0, to) = CMCW - CO + ™), 

so that only the case m > n needs to be considered. In particular, for m = n we 
find 

(6.4) T {n,Q 1 n)= 1 -C 2 {n)- l -a2n). 

At even weight N , the previous identities leave exactly N/2 — 2 Tornheim sums 
of the type T(m, 0,n) unevaluated. In effect, the convergence conditions require 
n > 2, so that all the sums with to = N/2 + 1, N/2 + 2, . . . , N— 2 are undetermined. 
The remaining N/2 — 2 unevaluated Tornheim sums are not all independent, since 
they satisfy linear relations obtained by applying Huard's identity (|1.4p to some 
special cases with known evaluations as, for instance, 

(6.5) T(m,k,0) = ((m)((k), 

{m 
T£(-iyt(i)((N-i)+ 

(6-6) £ C(i ) C (AT-i) (JV + 1)C(JV)^, 

i=2 



1 / N ~ 2 

(6.7) T(to, 1, 1) = - (JV + l)C(N) - V C(i)C(N 



\ i=2 / 

These results appear in (TTj . As usual, N is the weight of the Tornheim sum in the 
left hand side. 

On the other hand, Granville [9] showed that the multiple zeta values defined in 
satisfy 

(6-8) £C(Pi,P2,---,P fl ) = C(JV), 

where the sum is over all elements of Z ff such that pi + p-2 + ■ ■ ■ + Pg = N, with 
> 1 and pi > 2. In particular, when 5 — 2 we obtain 

(6.9) Yl T(m,0,n) = C(N), 

m+n—N 

where the sum is over pairs (to, n) with to > 1 and n > 2. 

Experimental observation. Not all of the relations stated above are linearly 
independent. However, it is possible to show that all Tornheim sums in for 
weights N = 4 and N = 6 can be completely evaluated in terms of zeta values (the 
details are provided at the end of this section) . In the case of even weight N > 8 
the Tornheim sums can be expressed in terms of zeta values and a reduced number 
of basis sums, which we choose as T(N — 2k, 0, 2k), with k — 1, . . . , K, where K is 
given by 

N -2 
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In particular, we need only one basis sum for weights 8, 10 and 12; two for weights 
14, 16 and 18, and three for weights 20, 22 and 24. The reader will find in [1] that 
K given above is an upper bound for the number of Tornheim sums required. 

Example 6.1. Indeed, for weight N = 8 we find that all sums in Zg have either 
explicit evaluations or can be expressed in terms of just T(6,0,2): 

T(0,0,8) = C(7)-C(8), 

r(i,o,7) = fc(8)-C(3)C(5), 

T(2,0,6) = §C(8)-T(6,0,2), 

T(3, 0, 5) = -^C(8) + 5C(3)C(5) + f T(6, 0, 2), 

r(4,0,4) = ^C(8), 

T(5, 0, 3) = 4f C(8) - 4C(3)C(5) - § T(6, 0, 2). 
Note that the identity (| 1 . 14[) follows directly from these expressions. 



Example 6.2. For weight N = 14 we can express all Tornheim sums as zeta values 
plus the sums T(12, 0, 2) and T(10, 0, 4): 

T(0,0,14) = C(13)-C(14), 

T(1,0, 13) - iiC(14) - C(3)C(H) - C(5)C(9) - W)\ 
T(2,0,12) = f|iC(14)-T(12,0,2), 

T(3,0, 11) = -3jjgk(14) + 11C(3)C(H) + 16C(5)C(9) + 9C(7) 2 + f T(12,0,2), 
T(4,0,10) = ^C(14)-T(10,0,4), 
T(5,0,9) = ^C(14) - ±f C(3)C(H) - 147C(5)C(9) - ^C(7) 2 

+ |T(10,0,4)- ±f T(12,0,2), 

T(6, 0, 8) = -^C(14) + 132C(3)C(H) + 240C(5)C(9) + 141C(7) 2 

-6T(10,0,4) + 66T(12,0,2), 

T(7,0,7) = iC(7) 2 -iC(14), 

T(8,0,6) = ^C(14) - 132C(3)C(H) - 240C(5)C(9) - 141C(7) 2 

+ 6T(10,0,4)-66T(12,0,2), 

T(9,0,5) = -«C(14) + ^ C(3)C(H) + 148C(5)C(9) + 2f C(7) 2 

-|T(10,0,4) + ifT(12,0,2), 
T(ll,0,3) - ^C(14) - 10C(3)C(H) - 16C(5)C(9) - 9C(7) 2 

-£T(12,0,2). 

We are now in a position to formulate a systematic and exhaustive algorithm 
to evaluate all the Tornheim sums in Zn. The reader is invited to download 
the Tornheim Mathematica 6.0 package developed by the authors, available at 
http://www.math.tulaiie.edu/~vnm/packages.html. Most of the calculations in 
this paper can be easily reproduced with the aid of this package. 
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Algorithm. The process of determining the Tornheim sums T(m,k,n) proceeds 
as follows. 

Step 1. First catalogue all sums that have known explicit evaluations. All of these 
cases appear in In this category we find: 

The sums with third entry n = 0: 

oo oo 

(6.11) T KM) = EEl^fe =C(m)C(fc). 

r— 1 s— 1 

The sums with third entry n = 1: 

f m 

T(m, fc, 1) = (-1)" 1 ^(-lJ'CCOC^ - *)+ 



, i=2 



(6.12) 

V-2 

-2C(i)C(iyr-i)--(JV + i)C(A0^ 



i=2 

where A = to + k + 1 is the weight. 



The sums of the type 

n 

(6.13) T(l, 1, n) = (n + l)£(n + 2) - £ C(0C(" + 2 - *). 



i=2 



The following sums in Zjy: 
(6.14) T(0,0,n) = C(n-l)-C(n), 

(6.15) 



T(l, 0, n) = i ^<(n + 1) - ^ C«C(n + !-*)! 



and the symmetric sum, 

(6.16) T(m, 0, to) = ^C 2 M - -C(2m), m > 2. 

For weight N = m + n odd, the sum T(m, 0, n) is given by 

TO — 1 



(6.17) r(m 1 0,n) = (-ir^ ( N ^ 1 )c(2j)C(N - 2j) 



Lf-J 



(-!) m E ( 7V n 2j i 1 )c(2i)c(^-2i)-KW- 



Step 2. If fc 7^ 0, use the reduction of Huard et al. given in (|1.4| to reduce the 
T(to, fc, n) not covered by the previous step to a finite sum of Tornheim sums in : 
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(6.18) 

~, , n v^/f+l;-!-l\ m , , v-—v (m + k — i — 1\ „. 

T(m,k,n) = Y / [ )T(i,0,N-i) + Y / [ ._. )T(i,0,N-i), 

i=i v m * ' i=l v 1 / 

with AT = m + k + n. 

Step 3. For N = m + n < 6 even and m, n > 2, compute all the sums T(m, 0, n) 
explicitly by solving the set of simultaneous equations obtained from: (a) applying 
Huard's reduction of the previous step to the identity (]6 . . for all independent 
pairs (m, k) with m+k = N; (b) Euler's identity (|6.3j) : and using the known explicit 
evaluations of the sums in Z 1 ^ given in Step 1. 

Step 4. For N = m + n > 8 even and m, n > 2, we write all the sums T(m, 0, n) 
in terms of the irreducible basis for weight N, 

(6.19) ^T(N-2k,0,2k), k=l,..., 

Step 5. The irreducible sums of the previous step are evaluated in terms of the 
integrals X^j using Theorem 12.31 

For a given even weight N = m + k + n, the whole process gives T(m,k,n) 
as a finite sum (with rational coefficients) of the zeta values ((N) and ((N — 1), 
products of two zeta values of the form C(j)((N — j) with 2 < j < N — 2 and, for 
N > 8, a finite number of integrals the type Y£ r N _ 2r , where 

(6.20) Y* r<N _ 2r := Y 2r , N _ 2r + (-l)^" 1 P ~ ((N - 1) log27r. 

Definition 6.1. We say that (m,k,n) is an admissible triple if m > k and m,k,n 
satisfy the conditions rn + n > 2, k + n > 2 and N > 3 for the convergence of 
T(m, k, n). 

We present now the results for small weight N = m + k + n. The cases of weight 
3 and 4 are straightforward, as all admissible triples corresponds to cases where the 
Tornheim sum has an explicit formula. 

Weight 3 

The admissible triples are (0, 0, 3), (1, 0, 2) and (1, 1, 1). We obtain 

T(0,0,3) = C(2)-C(3), 
T(1,0,2) = C(3), 
T(1,1,1) = 2C(3), 

directly from (|6.ip . (|6.2p and (|6.13p . respectively. The reader will find in 3 many 
proofs of the last identity. 

Weight 4 

The six admissible triples are 

(0,0,4), (1,0,3), (1,1,2), (2,0,2), (2,1,1) and (2,2,0). 



N-2 
6 
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We obtain 

T(0,0,4) = C(3)-C(4)=C(3)-^ ) 

T(l, 0, 3) = \ (3C(4) - C(2) 2 ) = |C(4) = 

T(l,l,2) = 3C(4)-C(2) 2 -ic(4) = ^, 

T(2,0,2) = 1(C(2) 2 -C(4)) = ^C(4) = ^, 

T(2,l,l) = i(5C(4)-C(2) 2 )=^C(4) = ^, 

T(2,2,0) = C(2) 2 -^C(4) = ^, 
directly from ([HE]), flO|), (|6~T3l> . (j^T^) . (j6~T2"j) and (f^TTT]) . respectively. 
The methods developed here produce the result 

(6.21) T(2, 0, 2) = ^ + 2 log(27r)C(3) - Y 2 , 2 . 
It follows that 

(6.22) y 2 , 2 = ^ + 21og(2^)C(3). 

As a consequence of this, we obtain the definite integral 
l 

(2ir 2 B 2 (q) Cli(27rg) - 27rBi(g) Cl 2 (27rg) - S (g) Cl 3 (2^)) lo g r(g) dq 

= ^+log(27r)C(3), 

where 

S (9) = l, B 1 (q)=q- 1/2, B 2 {q) = q 2 - q + 1/6, 

and 

m fo \ \^cos2irkq ^ sm2Trkq ^ cos2-nkq 

Cl 1 (2 7 rg) = ^ , Cl 2 (27T g ) = ^ , Cl 3 (27rg) = ^ — — . 

fc=i fc=i fc=i 

Considering that from (|5 . 32[) we also know that 

(6-23) J \ogT(q) Cl 3 (2^) dq = ~C(4) = 

we also find 

(6.24) f (2ir 2 B 2 {q)Ch{2irq)-2nB 1 {q)Ch{2irq))\ogr{q)dq 
Jo 



^+log(27r)C(3). 



Not your average integral. 



We state next the values of the Tornheim sums T(m, k, n) of even weight N > 
6. The formulas are the direct output of the algorithms presented here, the only 
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reductions used are those given at the beginning of this section. The reader will 
observe that the final expresions contain a single even value of the Riemann zeta 
function. This is artificial. For instance, the value 

T(4,2,0) = fC(6), 

should be written as 

T(4,2,0) = C(4)C(2), 

as in (|6. 1 1|) . This latter representation would be more helpful in the search for a 
closed-form expression for the Tornheim sums. However, at this point, we have 
decided to minimize the number of zeta values appearing in the formulas. 

Note 6.2. In the examples that follow, we do not list the Tornheim sum T(0, 0, ./V). 
This is the only sum that explicitly contains the term ((N — 1). 

Weight 6 

In this example we give complete details, which will be omitted for higher 
weights. The admissible triples are now (0,0,6), (1,0,5), (1,1,4), (2,0,4), (2,1,3), 
(2,2,2), (3,0,3), (3,1,2), (3,2,1), (3,3,0), (4,0,2), (4,1,1) and (4,2,0). 

A direct application of the identities and explicit formulas already discussed 
give the following evaluations (in the formulas below we have replaced the product 
C(2)C(4) by |C(6)): 

T(l,0,5) = |C(6) - K 2 (3), T(l,l,4) = |C(6) - C 2 (3), 

T(3,0,3) = -±C(6) + iC 2 (3), T(3,3,0) = C 2 (3), 

T(4, 1, 1) = |C(6) - iC 2 (3), T(4, 2, 0) = |£(6). 

Huard's expansion fll.4|) gives 

r(2,l,3) = 2T(l,0,5) + T(2,0,4), 

T(2, 2, 2) = 4T(1, 0, 5) + 2T(2, 0, 4), 

T(3, 1, 2) = 2T(1, 0, 5) + T(2, 0, 4) + T(3, 0, 3), 

T(3, 2, 1) = 6T(1, 0, 5) + 3T(2, 0, 4) + T(3, 0, 3). 

Euler's identity (jO)) . 

T(2, 0, 4) + T(4, 0, 2) = C(2)C(4) - C(6) = f C(6), 

allows us to express T(2, 0, 4) in terms of T(4, 0, 2). This, together with the explicit 
evaluations of T(l, 0, 5) and T(3, 0, 3), yields 



T(2,0,4) = 


|C(6) - 


T(4,0,2 


), 


T(2,l,3) = 


|C(6) - 


C 2 (3) - 


r(4,0,2), 


T(2,2,2) = 


|C(6) - 


2C 2 (3) - 


2T(4,0,2), 


T(3,l,2) = 


1C(6) - 


K 2 (3)- 


-2T(4,0,2), 


T(3,2,l) = 


7C(6) - 


|C 2 (3) - 


4T(4,0,2). 
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Finally, Huard's expansion (|1.4| applied to the Tornheim sum T(4, 2, 0) = |C(6) 
produces another identity, 

8T(1, 0, 5) + 4T(2, 0, 4) + 2T(3, 0, 3) + T(4, 0, 2) = |C(6), 

which permits to solve for the as yet undetermined value of T(4, 0, 2): 

(6.25) T(4,0,2) = f|C(6)-C 2 (3). 

This last result produces the explicit evaluation of all Tornheim sums of weight 6: 

T(l,0,5)= |C(6)-iC 2 (3), T(1,1,4) = |C(6)-C 2 (3), 

T(2, 0, 4) = -f C(6) + C 2 (3), T(2, 1, 3) = ±C(6), 

T(2, 2, 2) = ±C(6), T(3, 0, 3) = -|C(6) + ±C 2 (3), 

T(3, 1, 2) = -|C(6) + K 2 ( 3 )' T ( 3 > 2 ' !) = K 2 ( 3 )> 

T(3,3,0) = C 2 (3), T(4,0,2) = f|C(6)-C 2 (3), 

T(4, 1, 1) = |C(6) - K 2 (3)> T ( 4 > 2,0) = |C(6). 

All the Tornheim sums of weight 6 have been evaluated. 

Note 6.3. The problem of whether these sums are completely reduced is now equiv- 
alent to whether C 2 (3) and £(6) = 7r 6 /945 are rationally related. It is conjectured 
that £(3)/7T 3 is a transcendental numbei0. 



Weight 8 

The reduction algorithm described above begins by applying Huard's reduction 
procedure to express every sum T(m, k, n), with N = m + k + n = 8, in terms of 
the N — 2 = 6 sums 

(6.26) {T(l, 0, 7), T(2, 0, 6), T(3, 0, 5), T(4, 0, 4), T(5, 0, 3), T(6, 0, 2) } 
with k = 0. For example, 

(6.27) T(5, 2, 1) = 10T(1, 0, 7) + 5T(2, 0, 6) + 3T(3, 0, 5) + 2T(4, 0, 4) + T(5, 0, 3). 

We use the shorthand notation T(5, 2, 1) = [10, 5, 3, 2, 1, 0]. The table below gives 
all the coefficients corresponding to the 15 Tornheim sums of weight 8 with k ^ 0. 



1X1,1,6) = 


[2,0,0,0,0,0], 


T(2,l,5) 


= [2,1,0,0,0,0], 


T(2,2,4) = 


[4,2,0,0,0,0], 


T(3,l,4) 


-[2,1,1,0,0,0], 


r(3,2,3) = 


[6,3,1,0,0,0], 


T(3,3,2) 


= [12,6,2,0,0,0], 


T(4,l,3) = 


[2,1,1,1,0,0], 


T(4,2,2) 


= [8,4,2,1,0,0], 


T(4,3,l) = 


[20,10,4,1,0,0], 


T(4,4,0) 


= [40,20,8,2,0,0], 


T(5,l,2) = 


[2,1,1,1,1,0], 


T(5,2,l) 


= [10,5,3,2,1,0], 


T(5,3,0) = 


[30,15,7,3,1,0], 


T(6,l,l) 


= [2,1,1,1,1,1], 


r(6,2,0) = 


[12,6,4,3,2,1]. 







The authors wish to thank W. Zudilin for this information. 
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Therefore every Tornheim sum T(m, k, n) has been expressed in terms of the set 

(6.28) {T(i,0,N-i) : 1 < i < N- 2}. 

The value T(l,0,iV- 1) is given in and T(N/2,0,N/2) appears in (j6T)l . 

Moreover, Euler's relation (|6 . 3[) reduces the number of unknown Tornheim sums to 
N/2 - 2. In the case N = 8 the two unknowns are T(6, 0, 2) and T(5, 0, 3). Among 
the 15 sums discussed above, there are three with last index equal to 0, namely 
T(6, 2, 0), T(5, 3, 0) and T(4, 4, 0). Each one of them produces an equation in the 
unknowns T(6, 0, 2) and T(5, 0, 3) coming from the evaluation 

(6.29) T(m,fc,0) = C(m)C(fc)- 
For instance, the case T(4, 4, 0) gives 

(6.30) 5T(6, 0, 2) + 2T(5, 0, 3) = ^C(8) - 8C(3)C(5). 

This is the same relation among these sums obtained by Huard in Q1.14p . Unfortu- 
nately, the sums T(6,2,0) and T(5,3, 0) yield the same relation, so we are unable 
to produce an analytic expression for all Tornheim sums of weight 8, free of an 
unevaluated integral. 

We conclude that every Tornheim sum of weight 8 is a rational linear combination 
of the set 

(6.31) G 8 :={C(8), C(3)C(5), 7(6,0,2)}. 
The table shows the corresponding coefficients: 



T(l,0,7) 


= [f.-l.°] 


7(1,1,6) 


= [1,-2,0] 


T(2,0,6) 


= [5.0,-1] 


7(2,1,5) 


= [f,-2,-l 


T(2,2,4) 


= [f>-4,-2] 


7(3,0,5) 


- L 24 '°' 2l 


T(3,l,4) 


= [ — T'^' f] 


7(3,2,3) 


= [41 _! _1 

L 24 ' ' 2 


T(3,3,2) 




7(4,0,4) 


= [^,0,0] 


T(4,l,3) 


_ f 109 o 3] 
L 24 ' °' 2 J 


7(4,2,2) 


= [-f=2,l] 


T(4,3, 1) 


= [n.°.°] 


7(4,4,0) 


= [I,o,o] 


7(5,0,3) 


- [163 _ 4 „51 
L 24 ' ' 2J 


7(5,1,2) 


= [i-w] 


T(5,2, 1) 




7(5,3,0) 


= [0,1,0] 


7(6,0,2) 


= [0,0,1] 


7(6,1,1) 


= [I,-i,o] 


7(6,2,0) 


= [§.o,o]. 







The results derived in this paper give us the remaining unevaluated Tornheim 
sum T(6, 0, 2) in terms of the integral Y^q as 



(6.32) T(6, 0, 2) = |C(8) - 6C(3)((5) - K 



2,6- 
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Therefore, the generating set for Tornheim sums of weight 8 can also be taken 

as 

(6-33) G* s := {C(8), C(3)C(5), ^el- 



Weight 10 

For weight 10, the algorithm follows step by step the previous case. We find that 
all Tornheim sums in Z® are generated by the set 

(6.34) G w := {((10), ( 2 (5), C(3)C(7), T(8, 0, 2)}. 
For example, 

T(3, 2, 5) = -ifC(lO) + C 2 (5) + C(3)C(7) + ±T(8, 0, 2). 

According to Theorem 12. 3[ the remaining unevaluated Tornheim sum T(8, 0, 2) can 
be expressed in terms of the integral Y 2S as 

(6.35) T(8, 0, 2) = §C(10) - 8C(3)C(7) - 4C(5) 2 + Y 2 * s . 
In particular, we may also use the generating set 

(6-36) G* w := K(10), C 2 (5), C(3)C(7), Y 2 * s }. 

Generating set for Tornheim sums of even weight 

The same algorithm described above can be used to produce a generating set for 
the Tornheim sums of even weight N. Except for T(0,Q,JV) = ({N - 1) - £(JV), 
every such sum is a rational linear combination of the elements of the set 

(6-37) jc(iV), ((j)((N - j) : j odd , 3 < j < 2[^-\ + 1 

and, for N > 8, one must also include the collection of integrals 

(6-38) fe,iv-2r : 1 < r < X ~ 2 



where where Y 2 \ N-2r ' IS defined in (|1.22p . 



The smallest weight for which one requires two basis Tornheim sums is 14. In 
this case the generating set is 

(6.39) G u := {C(14), C 2 (7), C(5)C(9), C(3)C(H), T(12, 0, 2), T(10, 0, 4)}. 
The evaluation of T(12, 0, 2) and T(10, 0, 4) according to Theorem [231 gives 

T(12, 0, 2) = ||C(14) - 12C(3)C(H) - 12C(5)C(9) - 6C(7) 2 + Y* 12 , 
T(10, 0, 4) - £C(14) - 120C(3)C(H) - 60C(5)C(9) - 20C(7) 2 + F 4 * 10 , 

with 

^2*12= ^2,12 + 12C(13) fog 27T, 

yj; io = n,io+220C(13)log27r. 
In particular, we may also use the generating set 

(6.40) G* u := {((14), C 2 (7), C(5)C(9), C(3)C(H), Y 2 * 12 ,Y? 10 }. 



:S6 



OLIVIER ESPINOSA AND VICTOR H. MOLL 



Note 6.4. In [T] the reader will find the sums 

oo n— 1 

(6-41) ^,t):=EEw. 

n=l k=l 

that can be expressed as 

(6.42) a h (s, t) = C(s)C(t) - C(s + 1) - T(t, 0, a). 

The authors analyze a system of equations for the sums ah(s,t) with the weight 
w := s + t fixed. For w odd, the system has full rank and they obtain Huard's 
expression for the Tornheim sums. In the case w — 2n even, they establish that the 
dimension of the null space is \_(n — 1)/3J , so every sum can be expressed in terms 
of this number of basis elements. In our case, N = n/2, thus the expected number 
of generators for all Tornheim sums of weight N is at most [(N — 2)/6j . The fact 
that the set 

(6.43) { T(N - 2r, 0, 2r) : 1 < r < [(N - 2)/6j } 

can be used to generate all Tornheim sums (aside from the usual product of zeta 
values) will follow from a careful analysis of the identities generated by the relations 
T(m, n,0) = ((m)((n). We leave the details for the ambitious reader. The fact 
that these sums are linearly independent is beyond our reach. 



7. Conclusions 

We have discussed an algorithm that evaluates all the Tornheim sums 

oo oo _^ 

( 7J ) T(m,k,n) :=EE fms t ( r + sY^ 

r=l s=l V ' 

of a given even weight N := m + k + n, as rational linear combinations of the 
value C(-^0j products of zeta values ((j)C(N — j) with j odd in the range 3 < j < 
+ 1 and the integrals 



4 

(7.2) 
where 

(7.3) Y* 



N-2 
6 



2{2ir) m+n - 2 



m!(n-2)! ^ 



x 



j,m-\-n — 2—j 



N/2-1 [N-2 
m — 1 



C(iV-l)log27r, 



and the integral Xk t i is defined by 

/! 



(7.4) 



Xk,i := (-1) 



LZ/2J . 



(27T) 



\ogT (q)B k (q) Cl i+1 (27rg) dq. 



Here Bk is the Bernoulli polynomial and CI is the Clausen function. 



All the Tornheim sums of a given even weight N can be expressed in terms of 
zeta values and a reduced number of basis sums of the type T(N — 2r, 0, 2r), with 
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r = 1, . . . , [ ^g 2 j . These sums, in turn, can themselves be expressed in terms of 
zeta values and the integral i^* j^_ 2r , according to Theorem 12.31 




Our results may perhaps be used to develop fast numerical codes to compute 
even weight Tornheim sums to high accuracy. Since the whole family of Tornheim 
sums of a given weight can be expressed in terms of zeta values and a small basis of 
Tornheim sums, it is enough to compute the basis sums to the required accuracy. 
This will involve the numerical calculation of the Y integrals. 

For example, of the 46 admissible triples at weight N — 12, 30 give Tornheim 
sums that depend on the value of the single T(10, 0, 2). This sum is given by 
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(7.7) 



r(10, 0, 2) = 0.645 324 784 017 496 594 071 783 081 476 ... . 
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